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Abstract 

A linear spatial instability model for multiple spa- 
tially periodic supersonic rectangular jets is solved using 
Floquet-Bloch theory. It is assumed that in the region 
of interest a coherent wave can propagate'. For the case 
studied large spatial growth rates are found. This work is 
motivated by an increase in mixing found in experimental 
measurements of spatially periodic supersonic rectangu- 
lar jets with phase-locked screech and edge tone feedback 
locked subsonic jets. The results obtained in this paper 
suggests that phase-locked screech or edge tones may pro- 
duce correlated spatially periodic jet flow downstream of 
the nozzles which creates a large span wise multi-nozzle' 
region where a coherent wave can propagate. The large' 
spacial growth rates for eddies obtained by model cal- 
culation herein are related to the increased mixing since 
eddies are the primary mechanism that transfer energy 
from the mean flow to the large turbulent structures. 
Calculations of spacial growth rates will be presented for 
a set of relat ive Mach numbers and spacings for which ex- 
perimental measurements have been made. Calculations 
of spacial growth rates arc presented for relative Mach 
numbers from 1.25 to 1.75 with ratios of nozzle spacing 
to nozzle width ratios from s/icjv = 1 to s/ <e/v = 13.7. 
The model may be of significant scientific and engineering 
value in the quest, to understand and construct supersonic 
mixer-ejector nozzles which provide' increased mixing and 
reduced noise'. 


I. Introduction 

Interest in proving the economic and environmental 
feasibility of a high-speed civil transport has stimulated 
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studies of mixing enhancement in lobed mixer-ejector 
nozzles. By enhancing mixing the ejector length can be 
reduced with t he same' amount of noise suppression. In 
order to obtain information on such flows simpler configu- 
rations are studied. In particular, a simple' mixer nozzle 
configuration consisting of multiple rectangular nozzles 
with a synchronized screech instability was studied by 
Taghavi and Raman 1 and Raman and Taghavi 2 * 3 . This 
nozzle showed increased mixing with the jets synchro- 
nized. This paper use's the geometry and flow conditions 
investigated by Hainan and Taghavi 3 . The 1 same behavior 
is shown in a study of the effect of edge tones on multiple 
jet mixing of high-speed subsonic flows by Krothapalli et. 
al. 1 using the nozzle described by Krothapalli et. al. 5 . 

It is proposed that at some point before the jets merge 
local coherence can he achieved due to external forcing 
bv screech or acoustic feedback and large-scale propaga- 
tion of instabilities occurs with vary high growth rates. 
The temporal dynamics produced by the collective inter- 
action of compressible jets is discussed by Miles** and the 
collective' interaction of incompressible jets is discussed 
by Miles 7 . However, the predictions of spatial instabil- 
ity theory have shown better agreement with experiment 
in free shear flows and jets. Consequently, this paper 
presents a spatial instability analysis. Following Caster 8 , 
it is assumed the spatial instability analysis applies in 
a region when* nonlinear effects are small. However, it 
is acknowledged that the spatially growing linear theory 
will fail since the* amplitude of the disturbance must be 
bounded. 

For single nozzles a reduction in mixing and growth 
rates with increasing Mach number has been demon- 
strated experimentally by many investigators 9 13 . Cor- 
responding linear stability analysis of single nozzles 
shows results that are similar to the experimental 
studies 11 16 . This is attributed to the fact that eddies 
are the primary mechanism that transfer energy from 
the mean flow to the large turbulent structures. How- 
ever. the following study is bast'd on the idea that these 
experimental and theoretical results do not apply to the 
mixing of multiple supersonic rectangular jets with phase' 
locked screech. This paper is bast'd on a linear stability 
analysis of compressible periodic parallel jet flows which 
was undertaken to obtain results related to lobed mixer 
nozzles. In this study, the lobed nozzle design concept is 
extrapolated in a out' dimensional manner to arrive at an 
array of parallel rectangular nozzles separated by a dis- 
tance .s where t he smaller dimension of each nozzle is u?\ 
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and the longer dimension b is taken to be infinite. Note 
that it is assumed that even widely spaced rectangular 
jets which are phase-locked by screech are coherent spa- 
tially at some distance from the nozzle*. Consequently, in 
this linear stability analysis it is the collective behavior 
of compressible periodic parallel jet flow that determines 
the nozzle interaction. 

In this paper, the behavior of the solutions is discussed 
and the trace 1 of solutions is presented for a range of am- 
plification curves. For each operating condition, only the 
most highly amplified mode is of interest and special at- 
tention is paid to finding a good solution having the max- 
imum amplification. Calculations of spatial growth rates 
are presented for relative Mach numbers from 1.25 to 1.75 
with ratios of nozzle spacing to nozzle width ratios from 
t s/u’A r = 4 to s/icjsi = 13.7. The actual values are those 
for which experimental data is presented by Raman and 
Taghavi 3 . 

II. Results 

The nozzle configuration is shown in Fig. 1. In this 
paper, the flow is compressible and the velocity profile 
is adapted from an equation used by Monkewitz 1 ' in a 
study of the absolute and convective instability of two- 
dimensional wakes. A discussion of the problem formu- 
lation is given in Appendix A. Velocity profiles for ratios 
of nozzle spacing to nozzle 1 width of s/wjs; = 5.5 and 
& j w j\ r = 7.5 are shown in Figure 2. 

The linear spatial stability analysis is done using 
Floquct-Bloch theory. It is assumed that in the region 
of interest a coherent wave 1 can propagate. This type 
of analysis for temporal stability has been applied by 
Beaumont 18 to an incompressible flow with a sinusoidal 
velocity profile perpendicular to the flow and by Miles 6 *' 
to a compressible arid incompressible periodic parallel jet 
flow. This analysis procedure is discussed in Appendix 
B. 

Stability information is obtained using the flow model 
described in Appendix A and the Floquet-Bloch method 
described in Appendix B. In this study of spatially grow- 
ing waves proportional to exp {t(kx — u)r)} where k = 
kL*. x — x/L + . Cj — ljL*/AU and r = lAU/L * ,the flow 
disturbance is characterized by a real frequency, Cj. and 
a complex relative phase velocity , c = c r + %c t where 
Cj/k ~ c/AU = U/AU 4- c/2. Consequently, the phase 
velocity eigenfunction, c r , represents the phase velocity 
scaled by AU/2 and shifted by U where AU = f/ 2 — U\ 
and U = (U[ 4- f/ 2 )/2 so that for c r = 1 the disturbance 
moves at velocity f/ 2 , for ty = —1 the disturbance moves 
at velocity Ui and for c r = 0 the disturbance moves at 
velocity (/. For a given value of jet, relative Mach num- 
ber, m 2 , a value of the ratio of mean velocity to velocity 
difference. U/AU. a ratio of inter jet spacing to rectan- 
gular nozzle smallest dimension, {s/wj y). and c t . a range 1 
of Cj are studied to determine if a growing disturbance 


with amplification ,-A t , characterizes! by a periodicity 
parameter T t and a convective phase velocity c v exists. 
The coini>uter program first evaluates solutions at one 
hundred fixed values of <y in the range — 1 < c r < 1. The 1 
wave number is given by k — Cj/ {U/AU 4 c r / 2 4 A/2). 
Positive values of c t are used in the range from 0 to 1 at 
intervals of 0.1, since this produces negative k 7 values in 
a useful range. A solution at a given value ry is tabulated 
if the calculated value of F r is smaller than 5 .E — 04 and 
the calculated value of \S\ ( defined in Appendix B) is less 
than 2. A further search is made in the ry region where 
T r is smallest to find the best value of ry. An acceptable 
solution has 1 7 smaller than l.E — 06 and the calculated 
value of \6\ less than 2. 

The reported results at each value of Cj are limited to 
three : no solution, one solution, or two solutions. It is 
possible that more than two solutions exist. The model 
was developed to study amplification . —A,, over a range 
of Mach numbers and flow geometries for compressible 
periodic parallel jet flow when the flow is correlated be- 
tween the jets. 

The stability model is for shock-free supersonic jets 
where no screech tone exists. However, it does depend 
on the presence of a large span wise multi-nozzle region 
where* a coherent wave can propagate. In this paper, 
it is suggested that this region can be created by phase 
locked screech or edge-tones. Since screech generally oc- 
curs within a frequency range where the instability waves 
are highly amplified, the results from this study are used 
to explain certain events in screech synchronized multiple 
jets. 

For each condition studied, solutions for a range of Cj 
values at a given value of Cj were produced to find the 
region where the maximum amplification. — k 7 , of the un- 
stable wave* occurred. The value of r ? used were between 
0. and 1. using steps of 0.1. The value of Cj used was ini- 
tialized at -0.005 and incremented by 0.005. In general, 
blocks of 50 Cj points were examined at. a one time and 
the calculation for a particular value of r, was abandoned 
if the current block of 50 points and the previous block 
of 50 points had no solutions. 

To provide 1 information on the spatial instability solu- 
tion space, the trace of solutions for calculations of spa- 
tial amplification is presented for U/AU — 0.5. and a 
range of relative Mach numbers. m 2 = AU /a?* and of 
and nozzle spacings s/w n shown in Tables 1 and 2 where 
the values selected are those for which experimental data 
is presented by Raman and Taghavi 3 . 

Figures 3 thru 8 show plots of phase velocity eigen- 
value, ry , amplification, k t . periodicity factor. T, as a 
function of frequency, Cj for 0.1 < <: t < 1 using steps of 
0.1 for c,. For a given value of c x the trace of points for the 
spatial stability solutions can be characterized as having 
two regions. At low frequencies the trace shows region or 
band where the unstable solutions are continuous. This 
region is generally followed at higher frequencies by iso- 
lated islands of instability. 
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Parameter values for good solutions at the maximum 
growth rate for each case are presented in Tables 1 and 2. 
Values are selected so that T 7 . is smaller t han 1 .E - 06. 
Some points in Figure 71) and 8b have larger values of 
—kj. However, for these points the value of Y r was larger 
than \.E — 06 and they were 1 rejected. In order to find a 
good solution a progression of points was examined until 
an acceptable point with F r is smaller than 1 .E — 06 was 
found. 

The nozzle width, ws*, is 0.0069m. The frequency of 
the instability is given by 


uj r m-2(i o 


whore «o I s the ambient speed of sound (nominally 333 
m/sec). Values are given in Tables 1 and 2. Also shown 
in Table’s 1 and 2 is the excitation screech frequency. 
fg. The predicted instability frequency is about half the 
excit at ion frequency. 


III. Discussion 

A summary of the results of a study of temporal growth 
rates by Miles 6 is shown in Figure 9. Large* growt h rate’s 
were found for a range of spacings ,.s /)/,>'. and relative 
Mach numbers, m *- At larger spacings and higher Mach 
numbers the temporal growth rate is reduced. 

Figure 1 10 shows spatial growth rates. This paper 
shows similar trends when the spatial growth rate* is cal- 
culated. 


IV. Concluding Remarks 

A linear instability model for a large span wise multi- 
nozzle region far downstream where a coherent- wave can 
propagate is presented. Multiple supersonic rectangular 
jets exhibiting phase-locked screech or excited by edge 
tones may create such a region. The model may explain 
an increase in mixing observed in multiple jets exhibit- 
ing phase locked screech. This work was conducted with 
the expectation that multi-jets with synchronized screech 
could provide increased mixing and reduced noise. 

It might be that phase locked screech or edge tones 
can provide* a confining mechanism which produces spa- 
tial coherence* just as neighboring jets provided a con- 
fining mechanism in the experiments of Villermaux and 
Hopfinger 19 and VilleTinaux, Clagnc, and Hopfinger*“ u . 

The* model may be of significant scientific and emgi- 
neering value* in the quest to understand anel construct 
supersonic mixer-e v joctor nozzles. 


Appendix A: Formulation of the problem 

Let (U(y)A)A)) be the velocity of a steady plane- 
parallel How, where the x-axis is in the direction of the 
flow anel 


i'(y) = V • 4 ^- /((</) 

where* U\ is the velocity outside the jet. U> is the mean 
centerline jet velocity, U = 1 1 \ 1 2 . t he velocity scale* is 
A(/ = U -2 — U j, and /;.(//) is the velocity profile function 
which varies from - I to 1. 

The flow fie*ld is perturbed by introducing wave dis- 
turbance’s in the velocity and pressure wit h amplitudes 
that are a function of y. These disturbances are assumed 
to be traveling wave’s that are* expre*sse*d in dime*nsionle*ss 
variables as 


(w,f\ u\p) 

= ( w(j? ). HU)- p [* (o - . 

where* f/, f\ and tr are* dimensionless velocities, p is the* 
pressure*, r is tlie* time*, x is the distance along the jet 
flow, and y is the distance normal to the je*t along the* 
row of nozzles. The dimensionless variable* use*d herein 
are x = x/L *. y — y/L *. and r = tAU/L* 

The following quantities are also non-dimensionalizcd by 
a length scale, L* and a velocity scale, AU. 


A- = kL\ 

u JL* 


A U ‘ 

ijJ uJ 

T- = k±U 

and we define* r as follows 

" ir 


r 

AU 


U 


AU ^ 2 


By definition fex spacial instability, the frequency, A. 
is a real positive* number, k is a complex number that 
represents the wavenumber in the* x-dire*ction. —k t is the 
amplification rate of the* disturbance in the x direction, 
c a complex number that re*pre*se’nts the relative* phase* 
velocity. Since* 


by keeping c, positive* and varying e> betwee*n 4 1 and — 1 
ki is always negative* and only spatially growing solutions 
will lx* found. 

From the equations of motion if nonlinear and viscous 
terms are* neglected one* can obtain an equation for the 
y ~ component of the perturbation velocity as follows: 
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Appendix B: Floquet-Blocli theory 


f .4' 


- V 


h" x ? A* ti 

— 4- Aik — ( ~ -4 r ) — ~ 

-r) T A ( h-r ) 


= 0 ( 1 ) 


(h-r) T 

where the primes denote differentiation with respect to 

ti 

1 f 2 ; (h ” r ') 2 

A = —ik + m tk- 


4 

,4' = 2 


nr = 


f 


and from Crocco's Equation 21 

. T(y) T> (1 4 h(y)) 

T(y) = = — + 


Ti T v 


(i-S) 

^ I 


1, ,,, -f l)(/*(v)-l) 

~ 2 ( m i ~ O 4 


where 


MI 


m> = 


MI MJ y/T 2 y/Tk 

nil — “ n=p~ — pfr * 

«1 «2 v7i ViA 


where* a i is the local velocity of sound outside* the jet,, a *> 
is the local velocity of sound inside the jet. T v is the local 
temperature outside the jet, and T 2 is the local temper- 
ature inside the jet. 

In this paper, the velocity profile function. h(y). is pe- 
riodic such that 


h(y 4- 2ir) = h(y). 

The velocity profile h(y) is not any exact solution of the 
N a vier- Stokes equation, but it can be considered as a 
simple model of some real periodic flow. 

The velocity profile h{y) discussed herein is given by 

h(y) = 1 - 2 f(y) 


Since the basic flow velocity profile , f(y), is peri- 
odic. equation (1) is an example of a Floquet-Blocli 
problem. The mathematics of solving Floquet-Blocli 
type problems is discussed by Ince 22 . Hochstadt 23 , arid 
Zwillinger 21 . Applications to solid state* physics are dis- 
cussed by Sachs 25 . Brillouin 26 , and Dokker 2 '. Applica- 
tions to spatially periodic flow are discussed by Lorenz 28 , 
Green 29 , Beaumont 18 , Bai 30 and Gotoh 31 * 32 . 

The paper by Beaumont. 18 and the description of the 
Floquet-Blocli theorem by Hochstadt 23 were particularly 
useful in guiding this research. 

A survey of t he spatially periodic flow literature is pre- 
sented by Iv. Col oh and M.Y. Yamada 33 . 

The second order differential equation can be described 
by a system of first order different ial equations. Let 

v = Xj 


so that Eq. 1 can be r<‘writt.en as the system 


x, - jt > 


r', = Ds i -i Cx> 


( 2 ) 


where 


V A' , 

<=, r + X 


and 


D = 


h" 


(/» - r) 


Ad- -it A '~ h ' 


f .4 ' (h - c) 


If 4?(/y) is a fundamental matrix solution of equation ( 
(2 ) ) such that 


*(0) = I 

when* I is the identity matrix, then from the Floquet- 
Blocli theorem 


$(y + 2tt) = 4(y)*(27r) 


where the function f(y) is given by 

1 


f(y) = T- 


1 + ( siuh (^lf(T))) 


IS 


;/ = A ( — 1 4 £). and v goes from 0.0 to 2 tt. The 
profile function f(v) is adapted from an equation used 
by Monkewitz 17 in a study of the absolute and con- 
vective instability of two-dimensional wakes. Only two- 
dimensional disturbances will be considered. A schematic 
of the nozzle geometry is shown in Figure 1. Velocity 
profile using A = 1.22414 and A = 1.73897 are shown in 
Figure 2. 


We now introduce two solutions of equation ( 2 ) with 
initial values at y = 0.0. We have 4>(0) — [0i (0), <M0)] 
where <j>(( J)i — 1, 0'(O)i = 0, 0(O)*2 = d, and 0(O).> = 1. 
Next we seek the eigenvalues of #(2 tt) 

|*( 2 tt)-mI| 

= ft 2 — (0i(2/r) 4- ^>2 , (2n))/i 
+ (0,(27r)0 2 '(27r)- 
<h( 27T)0/(27T)) 

= p 2 ~ (0i (2 tt) 4 0 l >'(27t))// 4 1 = 0 (3) 
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Since 

^ 1 (27r)^ / (27r) - (2tt) = |4>(2tt)| - |*(0)| = 1 

The independent solutions of equation (2) have the 
form 

<p = X(y) cxp( ^|^ y) = A'(t/)exp{rj/) 

2tt 

The parameter T specifies the period of the eigenfunc- 
tion 4>. If T is real the eigenfunct ion grows or decays at 
infinity. Consequently, only imaginary values of I are 
acceptable. Thus the eigenfunction oscillates in space' 
and is called a continuous mode. The disturbance wit h 
T, = 1 /n. where n is a nonzero integer. has a period 2 / / tt . 
One with F, = 0 has the same period 2 tt as the main flow, 
while' an irrational value of F f means the disturbance is 
aperiodic. Note that the parameter I does not appear 
in the flow equation, but is due to the Floquet-Bloch 
theorem. 

Solutions of 2 are thus of the form 


X[(y - 2 tt ) = fi[Xi(y) 

X 2 (y - 2 tt ) = fi>X 2 (y) 

where pi and //•_> rc'present the zeros of (3). provided they 
are distinct. 

In general, these solutions will not be periodic. 

Conditions for periodic solutions can be found as fol- 
lows 

Let //j — ( 1( *‘ and /i L > — ( ~ }0i . 

Then from equation (3) 

cos(0/) = 4>\(2~) -*■ (p 2 (2n) — 6/2 


Consequent ly, for a solution to be periodic 6 must be real 
and \8\ smaller than 2. 

The constants // are t.ermed the characteristic multipli- 
ers of the Floquet-Bloch system (2) and the correspond- 
ing characteristic exponents are determined by the rela- 
tion T = I r 4 ilj = — + t^±. 
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TABLE I . — PARAMETER VALUES AT MAXIMUM 
GROWTH RATE FOR m 2 = 1.25, = 1.35 

AND m 2 = 1.45 


Parameter 

nij = 1.25 

m, = 1.35 

mj = 1 .45 | 

A 

1.294735 



shv„ 

4 

5.5 

11 1 

L *, meters 

gana* 



/,. Hz 

12192 

10200 


/*, Hz 

6696 

5412 

4651 

/,*//. 

0.549 

0.545 

0.546 

w r 

0.555 

0.54 

0.565 


-0.29605 

-0.29936 

-2.2653 


0.4 

0.4 

0.4 

k r 

1.192 

1.1625 

1.1527 

h 

-0.6773 

-0.6637 


6 r 

0.50985 

-0.90483 

-1.8499 

8 t 



1.4473x10 7 

r, 

0.2897 

0.17528 

6.2044xl0^ 2 

r. 




m 

0.50985 

0.90483 

1.8499 j 


TABLE II — PARAMETER VALUES AT MAXIMUM 
GROWTH RATE FOR m 2 = 1.55, = 1.65 

AND m 2 = 1.75 


Parameter 

rr^ = 1.55 

mj = 1 .65 

rrv, = 1.75 

A 

1.139382 

1.1258595 

1.111845 


10 

11.5 

13.7 

L*, meters 


1.3727x10 2 

1.6143x10 2 

/,. Hz 

7136 

5952 

5280 

/*, Hz 

3876 

3344 

2183 

/,*'/. 

0.543 

0.562 

0.413 

w r 

0.57 

0.525 

0.38 


-0.105 


-0.81295 


0.4 

0.4 

1.0 

kr 

1.0617 

0.87758 

0.13735 

k 

-0.47449 

-0.33654 

-0.73431 

5 r 

-1.8094 

-1.5363 


5, 



1.1 x 10 7 

r 


0.11059 

0.24607 

r r 


HK££l£U!ll 



1.8094 

1.5363 

4.9838 xlO* 2 j 
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Figure 1. — Nozzle configuration. 
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Figure 2.— Typical velocity profile. (A = 1.173897; 
s/wfg = 7.5 and A = 1.22414; (s/w^ = 5.5). 
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Figure 3. — (a) Phase velocity eigenvalue, c r , verses frequency, w, for 0.1 
< cj < 1.0. (s/w n = 4.0; A = 1.294735; U/AU = 0.5; m 2 = 1.25; T 2 = 1). 
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Figure 3.— (b) Amplification, kj, verses frequency, w, for 0.1 < q < 1.0. 
(s/w n = 4.0; A = 1.294735; U/AU = 0.5; m 2 = 1.25; T 2 = 1). 
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Figure 3.— (c) Periodicity factor, Tj, verses frequency, w, for 0.1 < q < 
1.0. (s/w n = 4.0; A = 1.4735; U/AU = 0.5; m 2 = 1.25; T 2 = 1). 
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Figure 4.— (a) Phase velocity eigenvalue, c r , verses frequency, w, for 0.1 
< q < 0.9. (s/w n = 5.5; A = 1.22414; U/AU = 0.5; m 2 = 1.35; T 2 = 1). 
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Figure 4.— (b) Amplification, kj f verses frequency, w, for 0.1 < 6 \ < 0.9. 
(s/w n = 5.5; A = 1.22414; 0/AU = 0.5; m 2 = 1.35; T 2 = 1). 
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Figure 4.— (c) Periodicity factor, Tj, verses frequency, w, for 0.1 < q < 
0.9. (s/w n = 5.5; A = 1.22414; U/AU = 0.5; m 2 = 1.35; T 2 = 1). 
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Figure 5.— (a) Phase velocity eigenvalue, c r , verses frequency, w, for 0.1 
< Cj < 1.0. (s/w n = 7.5; A = 1.173897; U/AU = 0.5; m 2 = 1.45; T 2 = 1). 



Figure 5.~(b) Amplification, kj, verses frequency, w, for 0.1 < Cj < 
1.0. (s/w n = 7.5; A = 1.173897; U/AU = 0.5; m 2 = 1.45; T 2 = 1). 
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Figure 5.— (c) Periodicity factor, T\, verses frequency, w, for 0.1 < Cj < 
1.0. (s/w n = 7.5; A = 1.173897; 0/AU = 0.5; m 2 = 1.45; T 2 = 1). 
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Figure 6— (c) Periodicity factor, Tj, verses frequency, w, for 0.1 < q < 
1.0. (s/w n = 10.0; A = 1.139382; 0/AU = 0.5; m 2 = 1.55; T 2 = 1). 
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Figure 7 — (a) Phase velocity eigenvalue, c r , verses frequency, w, for 0.1 
< Cj < 1.0. (s/w n = 11.5; A = 1.1258595; U/AU = 0.5; m 2 = 1.65; T 2 = 1). 
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Figure 7.— (b) Amplification, kj, verses frequency, w, for 0.1 < q < 
1.0. (s/w n - 11.5; A « 1.1258595; U/AU = 0.5; m 2 = 1.65; T 2 = 1). 
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Figure 7.— (c) Periodicity factor, Tj, verses frequency, w, for 0.1 < q < 
1.0. (s/w n = 11.5; A = 1,1258595; U/AU = 0.5; m 2 = 1.65; T 2 = 1). 
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Figure 8.— (a) Phase velocity eigenvalue, Cf, verses frequency, w, for 0.1 
< q< 1.0. (s/w n = 11.5; A = 1.111845; Cl/AU = 0.5; m 2 = 1.75; T 2 = 1). 
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Figure 8.— (b) Amplification, kj, verses frequency, w, for 0.1 < Cj < 
1.0. (s/w n = 11.5; A « 1.111845; Cl/AU = 0.5; m 2 = 1.75; T 2 = 1). 
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Figure 8.— (c) Periodicity factor, T\, verses frequency, w, for 0.1 < Cj < 
1.0. (s/w n = 13.7; A = 1.111845; 0/AU = 0.5; m 2 = 1.75; T 2 = 1). 
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Figure 9.— Temporal growth rate, (coj)max' results 6 for 
a range of relative Mach numbers, m 2 , and nozzle 
spacings, s/w n . 



Figure 10.— Spacial growth rate, (wj)ma* results for 
a range of relative Mach numbers, m* and nozzle 
spacings, s/w n . 
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